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6 Abstract 

7 We used coherent light scattering in a multi-speckle detection scheme to investigate the mesoscale 

8 dynamics in aqueous foam. Time-resolved correlation of the scattered speckle intensities reveals 
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9 the details of foam dynamics during aging. We introduce Temporal Contrast Analysis, a novel 

10 statistical tool that can be effective in characterizing structural rearrangements. Using Tempo- 

11 ral Contrast Analysis we were able to detect two distinct dynamical components present during 

12 foam aging: spontaneous and intermittent, avalanche-like events and continuous, flow-like rear- 
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13 rangements in the foam structure. We were able to measure these contributions separately from 

14 the intrinsic statistical noise contribution, and thereby independently analyze the decay of each 
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15 dynamical component during foam aging process. 
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i6 I. DYNAMICS OF JAMMED SYSTEMS 



17 The dynamics of densely packed grains, particles, droplets and bubbles are of great in- 

18 terest in the field of soft condensed matter, with examples ranging from granular piles and 

19 colloidal dispersions to foams, emulsions and biological tissues. As the density of particles 

20 is increased, these systems are often found to exhibit kinetic arrest. This kinetic arrest is 

21 characteristic of a "jammed" state in which the motion of any given particle is restricted due 

22 to the large number, and close proximity, of neighboring particles. As the system approaches 

23 a jammed transition, particle motion becomes more collective in nature, since in order to 

24 displace a single particle many other particles must move as well. As a result, many jammed 

25 systems exhibit slow, often non-equilibrium, relaxation commonly characterized by dynam- 

26 ical heterogeneities. This relaxation process is punctuated by intermittent, avalanche- like 

27 structural rearrangements 111, in contrast to continuous, flow-like (independent) motion of 

nnn 

28 particles in the low density regime. [2J |3| |4J Here we investigate the non-equilibrium coars- 

29 ening dynamics of aqueous foam using light scattering. Aqueous foam is one of the simplest 

30 systems exhibiting many characteristics of a jammed state. Despite being a combination of 

31 liquid and gas, foam retains many structural properties commonly found only in solids (for 

32 example, the ability to maintain its shape), as well as many properties of granular materials 

33 (such as flow under shear and jamming at low shear). Aqueous foam is also an example of 

34 a non-equilibrium system which, once formed, undergoes a continuous coarsening process. 

the rearrangement 
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35 During this process the average bubble size grows with aging time and 

36 rate of the tightly packed bubble structure decays over time. 

E 

37 One of the main difficulties in directly coupling to intermittent, avalanche- like structural 

38 rearrangements [12] is that they are both temporally and spatially rare and unpredictable. 

39 We studied the dynamics of aqueous foam using a variation of Diffusing Wave Spectroscopy 

40 (DWS) where one records the changes in the coherent diffraction pattern (or speckle) formed 

nn 

41 by laser light scattered off the internal foam structure [13( [1^]. The coherent speckle pattern 

42 effectively becomes a "fingerprint" of the specific configuration of scatterers (in this case the 

43 bubbles and the soap ridges forming the foam): a change in the configuration of scatterers 

44 will modify the pha ses and intensities of interfering waves and therefore result in a different 

45 speckle pattern [l^] [3] [l^]. Therefore DWS is a highly suitable technique for studying foam 

46 coarsening since it provides a direct way of coupling to the dynamics of a large ensemble of 
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scatterers. 



48 II. MATERIALS AND METHODS 

49 The samples used here were Gillette Foamy Regular shaving foam. The typical dimension 

50 of the randomly distributed foam bubbles in this sample is on the order of a few micrometers, 

51 and the relatively long coarsening lifetime allows us to probe a substantial range of the 

52 dynamical life of the foam - from seconds to many hours. Due to the rapidly evolving 

53 nature of aging foam dynamics, we record man y sp eckle intensity values simultaneously, a 



54 so-called "multi-speckle" detection scheme 
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55 The coherent light source was a 2.0 mW HeNe laser, with a wavelength of 633 nm. The 

56 data were collected using a line scan charge-coupled device (CCD) detector. The 1024 pixel 

57 (pixel size 10 fim by 10 jttm) spL2048-140k Basler line-scan CCD camera was placed 6.2 cm 
5B from the sample at an azimuthal scattering angle of 45 degrees. We collected 32 batches 
59 of data at 1 minute intervals, starting 1 minute and 30 seconds after we first expressed the 
eo foam. Each batch of data consisted of 32,768 line-scan frames with 74 /is exposure time 
ei taken 1 ms apart, for a total duration of 32.768 seconds for each batch. 

62 The raw speckle data is shown in Fig. (TJintensity (indicated by color) over a 100 pixel 

63 CCD region is plotted as a function of time, over a total time period of 1 second, with four 

64 different panels corresponding to selected foam aging times of 1, 5, 15 and 30 minutes. For 

65 early aging times the speckle streak image appears chaotic due to the high number of rapid 
ee intensity fluctuations. As the foam ages the images show extended sections of constant 

67 speckle intensities, appearing as horizontal streaks in the image. These streaks represent 

68 the periods of time when the foam structure is fairly static. Periods of relative stability 

69 are occasionally interrupted by rapid rearrangements of speckle patterns; two examples of 

70 such intermittent rearrangements are identified with red arrows in the 30 minute aging time 

71 sequence (at approximately 240 and 550 ms). We first quantified the aqueous foam dynamics 

72 by calculating the intensity autocorrelation function as a function of foam age. Intensity 

73 autocorrelation function is defined as: 

92(r) = ^ (1) 

74 where I(t) and I(t+r) are speckle intensities recorded at time t and t+r, respectively, and 

75 averaging is performed over all accessible t and pixels while keeping r constant. Fig. [2] 
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FIG. 1. Speckle intensity (color scale) for 100 pixel CCD region (vertical axis) as a function of 
time (horizontal axis, increasing from left to right) over a one second time interval. Red arrows 
indicate sudden changes in speckle patterns. The characteristic (exponential) decay time constant 
increases with aging time. Each image has been scaled for visibility. 



76 shows the autocorrelation functions for 7 selected aging times (each was computed for a 32- 

77 second continuous batch dataset). As the sample ages, the characteristic decay time of the 
7B autocorrelation function increases, as expected. Fig. [3] shows the exponential relaxation time 
79 constants to obtained from the fits of the autocorrelation data shown in Fig. [2] to a fitting 
so model : ^(t) = a, ■ exp((r/toY). *, plotted as a function of aging time. Unlike the scaling 
si behavior for very early ages (less than 1 minute), when the rate of bubble rearrangements 
82 during coarsening is found to scale as t~ 2 Q , we observe a t ~ t~ °' 34 scaling over the range 
as of aging times from 1 to 500 minutes. 
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FIG. 2. Autocorrelation function at age 480 min. has been scaled by a factor of 5. 

85 III. TWO TIME CORRELATION 

86 A standard approach to characterizing non-equilibrium dynamics involves calculating the 

87 two-time correlation function 



24j (Eq. [2J, most conveniently displayed as a contour plot of 
83 92^1^2) a s a function of t\ and t%. As the dynamics of a system slow down (and characteristic 

89 relaxation time increases), the area of high correlation around the diagonal of the contour 

90 plot (where t\ = £2) will broaden with time. 



o(f . _ (1(h) ■ I(h)) q 

92(h,h)- (2) 

91 Based on the similarities between Eq. [1] and Eq. [21 one would normally expect the mono- 

92 tonic broadening of (72(^15 ^2) contours in an aqueous foam system as the foam ages. However, 

93 the two-time correlation function shown in Fig. H] exhibits a different behavior, with several 

94 regions of broadening with respect to ti = t 2 diagonal, separated by "pinch-off" points where 

95 the contour abruptly narrows. This behavior is caused by a substantial intermittent dynam- 

96 ical component in foam samples: whenever there is a spontaneous (or intermittent) global 
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FIG. 3. Relaxation times to obtained from fits of autocorrelation functions shown in Fig. [2l as a 
function of foam age. The best fit to the data is to = a • ta + b power law behavior, with 7 = 0.34, 
or roughly a cube-root power law. 



97 rearrangement in the foam structure, the entire speckle pattern changes and the correlation 

98 to earlier speckle patterns becomes, suddenly, very poor. Following the rapid rearrange- 

99 ment, the g^ifi: £2) contours broaden again due to continuous dynamics, until the next large 

100 intermittent rearrangement, at which point the contour of ^2(^1^2) pinches off again. Al- 
ien though two-time correlation analysis fails to illustrate the diverging characteristic timescale 

102 in this particular case, the shape of gv(t\, ^2) contours provides us with clear evidence of two 

103 distinct dynamic components of coarsening in aqueous foam, one continuous and another 

104 intermittent. 



6 




1 ?R 

1 . £.iJ 




1 7 


- 


1.15 




i i 
i . i 




1.05 






1 





50 100 150 
t (ms) 



200 



250 



FIG. 4. Two-time correlation function 52(^1^2) plotted for 250 ms time interval for 32 minute age. 



105 IV. TIME RESOLVED CORRELATION 



Another method for coupling to the intermittent component of dynamics in a jammed 
system is Time Resolved Correlation (TRC) analysis [1] . In contrast to Eq. (TJ which takes a 
time average of correlation values at a fixed time offset r, TRC (Eq. [3]) shows the correlation 
values for a single selected value of r as a function of time, t. TRC also corresponds to cross- 
sectional cut of the two-time autocorrelation function ^2(^1,^2) along a 45-degree diagonal, 
offset from the t\ = £2 line by ±r. 

(I(t)-I(t + r)) q 



TRC(t) 



(3) 



(I(t)) q -(I(t + r)) q 

For a dynamically homogeneous system the probability distribution of the TRC values is 
expected to follow a Gaussian distribution [1] . The TRC of the foam data at early aging times 
(less than 5 minutes) displays a Gaussian distribution. This is because foam rearrangements 
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FIG. 5. The probability distribution function of TRC values becomes asymmetric, skewed towards 
lower values as the foam ages, due to temporal inhomogeneities. 

us are happening so rapidly that the dynamics closely resemble that of a continuous system. As 

lie the foam ages, an asymmetric tail emerges at the lower end of the TRC distribution [Fig. |5]. 

H7 The distribution becomes skewed towards lower TRC values as the dynamics become more 

us temporally inhomogeneous. 



119 V. TEMPORAL CONTRAST 

120 While TRC provides clear evidence of heterogeneous dynamics with a substantial in- 

121 termittent component, the gradually varying baseline of TRC makes it difficult to clearly 

122 identify the times at which intermittent rearrangements occur. In order to better isolate the 

123 intermittent component we devised a new method called Temporal Contrast (TC). The con- 

124 trast (visibility) of spatially varying fringes is defined as the ratio of the variance of intensity 

125 with respect to the mean. We define a similar function as a ratio of variance to the mean 

126 for intensities measured within a specific pixel over a narrow time interval, r. We calculate 
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FIG. 6. Temporal Contrast (upper panel) and Time Resloved Correlation (lower panel) plotted 
over for 10 second period at a sample age of 31 minutes. Sudden structural rearrangements exhibit 
themselves as spikes in Temporal Contrast and coinciding dips in TRC data. 



the temporal contrast for each pixel, then average over all q (pixels). 

Imax Imin 



T.Contrast(t, r) 



(4) 



L max i ± min , 

Imax/min = max/min[I(t - r/2), ■ ■ • ,I(t), ■■■ ,I{t + r/2)] 

Where r describes the duration of time interval over which the contrast is calculated. In this 
study we used a r = 20 ms window which was selected as a compromise between statistics 
and the temporal resolution. By taking the average over q, continuous fluctuations are 
primarily suppressed and collective intermittent fluctuations show up as spontaneous spikes 
in Temporal Contrast. 

The example of Temporal Contrast shown in Fig. E] provides a much more detailed pic- 
ture of the sample dynamics - the spikes of various magnitudes correspond to intermittent 
structural rearrangements that vary in either extent of the spatial dimensions or the rela- 
tive displacements of the scatterers involved in the specific "avalanche event. As the foam 




FIG. 7. The frequency of intermittent rearrangements as a function of aging time. The fit to the 
data shown with a line is represented by an exponential decay with a time constant ij = 10.85 
minutes. 

138 coarsens the frequency of avalanche events clearly decreases and the Temporal Contrast 

139 value between the spikes stabilizes at a decreasingly lower, but finite baseline value. Fig. [7j) 
wo shows the frequency of intermittent, "avalanche-like rearrangements as a function of aging 
mi time. In calculating the avalanche frequency, we count the frequency at which Temporal 

142 Contrast value exceeded a certain threshold value, which was denned as the average be- 

143 tween the maximum contrast (mean of 50 highest spike values) and the overall mean of the 

144 Temporal Contrast in each 32 second exposure. 

145 One should note that if the sample structure is static for extended periods of time, the 
we scattering speckle pattern is also static, resulting in Temporal Contrast value that should 
147 approach zero. This is not the case for Temporal Contrast calculated for our data: for 
we example, in Fig. El the baseline value is about 0.04. There are two possible explanations 

149 for this finite value of Temporal Contrast: the first possibility is that finite photon counting 

150 statistics would lead to recorded intensity variations within each pixel even for absolutely 
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151 static sample, resulting in non-zero value of Temporal Contrast, while the second possibility 

152 is that the sample undergoes slowly varying, continuous structural changes which results in 

153 (relatively small) speckle intensity fluctuations. As we demonstrate below, both of these 

154 effects are present in our sample. When considering the effects of statistical noise from con- 

155 tinuous dynamics on the Temporal Contrast baseline value, one should note that statistical 

156 variations of intensities from their mean are generally completely random as a function of 

157 time, while intensity fluctuations due to continuous dynamics will follow a continuous and 

158 monotonic trend (decreasing or increasing at various rates for different pixels), and this 

159 rate of change , measured over a period of time much shorter than the characteristic re- 
leo laxation rate, will remain unchanged. In order to separate these two types of contributions 
lei we down-sample the raw intensity datasets by a factor n, producing a sequence of frames 

162 consisting only of every n-th exposure, and discarding the exposures in between. We then 

163 compute the Temporal Contrast baseline (defined as the mean of the 50 smallest Temporal 

164 Contrast values) in each of these down-sampled batches, and compare them as a function of 

165 down-sampling ratio n. 

lee If in between spontaneous "avalanche events the sample is truly static and the finite 

167 baseline value was due solely to the random noise present in the speckle, the down-sampling 

168 procedure will not change the baseline values of Temporal Contrast. However, if the sam- 

169 pies undergoes continuous structural evolution, resulting in continuous change in speckle 

170 intensities, the down-sampling will effectively compress these trends by a factor n, resulting 

171 in higher calculated baseline value of Temporal Contrast for these down-sampled datasets 

172 (note that Temporal Contrast will increase regardless of whether the intensity in the partic- 

173 ular pixel is increasing or decreasing). For small values of Temporal Contrast (the relative 

174 changes in pixel intensities are small) the rate of intensity changes due to sample dynamics 

175 can be assumed constant, which would result in the baseline value of Temporal Contrast for 

176 down-sampled datasets increasing linearly with n. 

177 Fig. [HJ shows the linear dependence of the Temporal Contrast baseline as a function 

178 of down-sampling ratio, n, for different aging times. By measuring the slope of the data 

179 with respect to n, we can directly measure the change in the continuous dynamics as the 

180 sample is aging as can be seen from Fig. [HJ the slope decreases during aging, which implies 

181 that the rate at which the sample undergoes continuous evolution is also decreasing over 
lea time. Figure [9] shows the baseline Temporal Contrast, as well as the rate of continuous 

11 




FIG. 8. The baseline value of Temporal Contrast calculated over each 32 second frame dataset 
for foam samples aged for 5, 8, 13 and 25 minutes as a function of down-sampling ratio, n, shown 
along with linear fits. The slope of the fit represents the contribution of continuous dynamics 
component, while the constant term (y-axis intercept at n=0) is due to the random shot-noise 
statistical contribution. 



184 dynamics level as a function of sample age. The difference between the baseline and the 

185 continuous component is the contribution to the baseline due to finite counting statistics, 
we or shot-noise. The continuous component of the dynamics decays exponentially over time, 
187 with a time constant comparable to the exponential decay time constant of the rate of 
las avalanche-like intermittent rearrangements, shown in Fig. [7J This remarkable similarity in 
189 behavior of intermittent and continuous dynamics may be indicative of a coupling between 
wo the spontaneous and rare avalanche-like (temporally heterogeneous) events and the slow 

191 continuous flow that occurs in between avalanches which corresponds to the temporally 

192 homogeneous component of the dynamics. One possible model of the coupling mechanism 

193 is that the flow-like motion occurs in response to stress fields and voids created as a result of 

194 sudden and drastic rearrangements of the foam structure that occur during and preceding an 
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FIG. 9. The change in baseline value of Temporal Contrast (red circles) and the rate of continuous 
dynamics (blue circles) obtained from fitting the slope of change in baseline for down-sampled data, 
as shown in Fig. [8] , as a function of aging time, t a . The difference between the two curves is the 
contribution to the baseline value of Temporal Contrast due to shot-noise fluctuations.. 



195 avalanche event. However, it is also possible that the flow-like continuous rearrangements 

we can themselves lead to building up of stresses that are ultimately culminated in a rapid 

197 avalanche event or series of avalanches, due to, for example, the bursting of soapy ridge or 

198 multiple ridges. 



199 VI. SUMMARY/CONCLUSIONS 

200 By using several traditional statistical methods of speckle intensity fluctuation analysis, 

201 we were able to show that dynamics in foam exhibit slowing behavior as the foam ages and 

202 coarsens. We also introduced a new method of studying these dynamics in foam, Temporal 

203 Contrast, which is especially effective in elucidating the intermittent, spontaneous structural 

204 rearrangements in systems consisting of a large number of scatterers. Through the statistical 

13 



205 down-sampling analysis of Temporal Contrast baseline values we were able to advance upon 

206 previous studies of foam dynamics by showing that in between intermittent rearrangement 

207 the foam sample undergoes continuous flow-like dynamics. The rate of this slow, continuous, 

208 dynamical component can be separated from the contribution of statistical noise to the 

209 Temporal Contrast data, and was shown to decay exponentially during sample aging, in 

210 striking similarity to exponential decay of the rate of intermittent rearrangements as a 
2n function of aging time. This may indicate that the continuous, flow-like rearrangements are 

212 coupled to the intermittent, avalanche- like events, even though the precise nature of coupling 

213 remains unknown. One can speculate, however, that the continuous flow-like dynamics may 

214 be due to the slowly varying internal stresses in response to foam drainage, while avalanche- 

215 like rearrangements are due to the bursts of the film connecting the bubbles. 

216 The statistical analysis methods of temporal speckle fluctuations applied by us to study 

217 foam coarsening will be useful for studying complex non-equilibrium dynamics in a wide 

218 variety of other systems, where both intermittent and continuous dynamics may be present. 

219 The examples of such systems include not only a wide range of soft matter and granular 

220 systems near jamming transitions, but also a wide range of hard condensed matter systems, 

221 such as dynamics of magnetic, ferroelectric, charge-ordered domains as well as systems in the 

222 vicinity of electronic, structural and magnetic phase transitions. The analysis of temporal 

223 fluctuations in visible light speckle is directly applicable to shorter length scales probed 

224 by techniques such as X-ray Photon Correlation Spectroscopy and Fluorescence Intensity 

225 Fluctuation Spectroscopy. 
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